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ABSTRACT 



Gamma-ray observations of a stellar tidal disruption event (TDE) detected by the Swift 
satellite and follow up observations in radio, mm, infrared and x-ray bands have provided a 
rich data set to study accretion onto massive blackholes, production of relativistic jets and their 
interaction with the surrounding medium. The radio and x-ray data for TDE Swift J1644H-57 
provide a conflicting picture regarding the energy in relativistic jet produced in this event: x- 
ray data suggest jet energy declining with time as t^^/'^ whereas the nearly flat lightcurves in 
radio and mm bands lasting for about 100 days have been interpreted as evidence for the total 
energy output increasing with time. We show in this work that flat lightcurves don't require 
addition of energy to decelerating external shock (which produced radio and mm emission via 
synchrotron process), instead the flat behavior is due to inverse-Compton cooling of electrons 
by x-ray photons streaming through the external shock; the higher x-ray flux at earlier times 
cools electrons more rapidly thereby reducing the emergent synchrotron flux, and this effect 
weakens as the x-ray flux declines with time. 

Key words: radiation mechanisms: non-thermal - methods: analytical - x-rays, radio: bursts, 
theory 
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1 INTRODUCTION 

When a star wanders too close to a massive blackhole it is shredded 
by the strong tidal gravity of the BH. In this process a fraction of the 
star is captured by the BH and is eventually accreted, and roughly 
half of the stellar mass is flung out on hyperbolic orbits (Lacy et 
al. 1982; Rees 1988; Evans & Kochanek 1989; Goodman & Lee 
1989; Ayal, Livio & Piran 2000). Such an encounter is expected to 
occur at a rate of 10~^ — 10~^ yr~^ per L» galaxy (Magorrian & 
Tremaine 1999; Wang & Merritt, 2004; Bower, 2011). The accre- 
tion disk is expected to produce blackbody radiation in UV and soft 
x-ray bands from the region close to the BH, and a bright optical 
flash is produced by a super-Eddington outflow and by the iiTa- 
diation and photo-ionization of unbound stellar debris (Strubbe & 
Quataert 2009, 2010). 

Generation of a relativistic jet is also expected to accompany 
a tidal disruption event (TDE) as the stellar material bound to the 
BH gravity is accreted over a period of time; we know empirically 
that accreting BHs in active galactic nuclei (AGNs) produce jets 
moving at speed close to that of light with Lorentz factor of order 
10. The jet launched in a TDE interacts with the circum-nuclear 
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medium (CNM) and produces synchrotron radiation over a broad 
frequency band from radio to infrared which Giannios & Metzger 
(201 1) predicted should be observable for several years (also Met- 
zger et al. 2012). 

The NASA Swift satellite discovered a transient event. Swift 
J 1644+57, in 7 & x-ray bands two years ago (Burrows et al. 201 1, 
Bloom et al. 201 1) which was subsequently observed in radio, mm, 
infrared and x-ray bands for a time period of more than year with 
good temporal coverage, e.g. Zauderer et al. (2011), Levan et al. 
(2011), Berger et al. (2012), Saxton et al. (2012), Zauderer et al. 
(2012), and the data is broadly in line with expectations of a TDE, 
e.g. Burrows et al. (2011). 

Two other transient events observed in the last year have also 
been suggested to be TDE (Cenko et al. 2011, Gezari et al. 2012). 
However, the follow up observations for these events have been 
sparse, and so we consider in this paper only Swift 1644+57 and 
the puzzle posed by the rich data set for this event. 

The mm & radio lightcurves for Sw J1644+57 were found to 
be flat/rising for ~ lO^days by Berger et al. (2012). These au- 
thors suggested that this implies continuous injection of energy 
to the decelerating blast wave, which produced this radiation, for 
~ lO^d. However, that requires about 20 times more energy in the 
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blast wavqj than the initial jet energy that produced the strong X- 
ray signal which was of order lO'^^erg. This energy requirement of 
more than a few times lO^^erg is very difficult to satisfy for a TDE 
model. We show in this paper that the multi-wavelength lightcurves 
of Sw J 1644+57 can be understood without requiring energy in the 
blast wave to increase with time. 

In the next section we describe the main features of the data for 
Swift 1644+57, and the main puzzles. We then describe in §3 the 
main idea of this paper that offers a simple solution to the puzzling 
data, and the application to Sw J1644+57 is provided in §4. Some 
concluding remarks about this event can be found in §5. 



& Kumar (2000), Granot and Sari (2002). Energy injection to the 
decelerating external blast wave between 10 and 10^ days by a fac- 
tor ~ 20 has been invoked, e.g. Berger et al. (2012), to prevent the 
decline of the radio lightcurve during this time interval. This huge 
amount of energy added to the blastwave, however, leaves no x-ray 
foot print, and that clearly is very puzzlingj. The model we suggest 
— IC cooling of electrons by x-ray photons — avoids this problem 
because it requires no energy addition to the decelerating external 
shock to flatten radio/mm lightcurves. 

Two other puzzles suggested by the data are the precipitous 
decrease of Va by a factor 10 between 5 and 10 days, and the rapid 
fall off of the infrared lightcurve during the same period. 



2 SWIFT J1644+57: SUMMARY OF OBSERVATIONS 
AND PUZZLES THEY HAVE POSED 

Swift/XRT found that the x-ray flux from Sw J1644+57 varied 
rapidly on time scale Stobs ~ 10^ s, indicating that the emission 
was probably produced at a distance of order 2c5toba^^j 5; lO^^cm; 
Vj < 10 is jet Lorentz factor. The average x-ray flux was roughly 
constant for the first 5 days, then dropped suddenly by a factor 50 
during the next 3 days, and subsequently for 20d < tabs H 500d it 

— 5/3 

declined as ~ t^^^^ (Burrows et al. 2011; Berger et al. 2012; Za- 
uderer et al. 2012). The x-ray flux then declined further by a factor 
~ 170 between 500 days and 600 days (Zauderer et al. 2012). The 
spectrum in the 0.3-10 keV Swift/XRT band is reported to be a 
powerlaw function — f^ (x u~^ — with index /? ~ 0.7 during the 
first ~ 50days and then decreasing to ~ 0.5 a few hundred days 
after the tidal disruption. 

This transient event received extensive coverage at multiple 
mm and radio frequencies. Depending on the frequency band the 
lightcurve was found to rise, or remained nearly flat, during the 
first 100 days, and then declined as roughly f^/^^. The spectrum 
at low frequencies ( < lO^^Hz) rose as ~ u^'^ which is roughly 
consistent with that expected of a synchrotron spectrum below the 
self-absorption frequency (Va)- And the spectrum above the peak of 
/^ fell off as ~ !/~'''^. The synchrotron self-absorption frequency 
(Va) underwent a dramatic decline by a factor of 10 between 5 and 
10 days, and subsequently Va displayed a very slow decline over 
the next several hundred days. The flux at the peak of the spectrum 
fell off during 5-10 days by a factor ~ 3, then increased by a factor 
^ 4 during 10-10^ days (Berger et al. 2012). It should be noted 
that these peculiar behavior of the radio data during the period 5- 
10 days coincided with a sharp decline of the x-ray flux by a factor 
~ 10^, e.g. Burrows et al. (2011), Berger et al. (2012). 

The event was also observed in the infrared K-band (~ 1.4 x 
lO^^Hz) multiple times between 5 and 50 days, and the flux was 
found to decline during the first 5 days as ^ t~^^ (Levan et al. 
2011). 

One of the main puzzles posed by the rich data for this event is 
that lightcurves in > lOGHz frequency bands, which lie above the 
synchrotron absorption and other characteristic frequencies, show 
a flat or a rising behavior when in fact according to the external for- 
ward shock model these lightcurves should be declining as ^ t~j,^ , 
e.g. Meszaros & Rees (1994), Chevalier and Li (2000), Panaitescu 



^ Using simple relativistic equipartition arguments Bamiol Duran & Piran 
(2013) have shown that in the context of a standard synchrotron scenario, 
without IC cooling of electrons, the requirement that the blast wave energy 
should increase with time during the first ~ 200 days in order to explain the 
radio lightcurve behavior, does not depend on the choice of the geometry of 
the jet, nor does it depend on the precise spectral fitting of the radio data. 



3 INFLUENCE OF X-RAYS ON EXTERNAL SHOCK 
RADIATION 

Basic idea: relativistic jet from the TDE interacts with the CNM 
(circum-nuclear medium) and drives a strong shock wave into ijj. 
CNM electrons accelerated by the shock radiate synchrotron pho- 
tons at radio and higher frequencies. These electrons are cooled by 
the IC scattering of x-ray photons streaming through the shocked 
plasmqj (these x-ray photons are the same radiation that was ob- 
served by the Swift satellite); see Fig. 1 for a schematic sketch 
of this scenario. Since the x-ray flux declines rapidly with time 
(fx oc t~^''^ for t > 10 days), the IC cooling of electrons also di- 
minishes with time and this is what is responsible for transforming 
an otherwise declining lightcurve (/^ Ztobs'^ ^° ^ ^^^ °^ ^ rising 
lightcurve in radio/mm band; the IC cooling by x-rays turns out to 
be very effective in suppressing synchrotron radiation in mm/cm 
bands at early times, and this suppression gets weaker with time as 
the x-ray flux declines. 

One other thing we should address here for those who might 



be thinking that the observed spectrum above the peak (/„ 



') 



is too shallow to be consistent with a fast electron cooling model. 
The observed parameters of Swift 1644+57 are such that IC scat- 
terings for electrons radiating above the peak takes place in the 
Klein-Nishina regime, and hence the shallow observed spectrum; 
detailed calculations are presented in §3.1 & 3.2. 

Since the main component of our proposed idea to explain 
radio/mm lightcurve behavior is IC cooling of electrons in the ex- 
ternal shock by the observed x-ray radiation, we begin the tech- 
nical part of this section with a description of this process. We 
show that it is unavoidable that IC cooling by x-ray photons signif- 
icantly modifies electron distribution in the external shock for Sw 
J1644+57. We also discuss whether IC scatterings take place in the 
Thomson or Klein-Nishina (K-N for short, from now on) regime. 



^ The x-ray lightcurve fell off as t ^ for t^bs ~ 5d, and therefore ac- 
cording to the x-ray data most of the energy release in relativistic jet took 
place during the first 5 days. 

•^ A reverse shock propagating into the jet can also be produced in the 
CNM-jet interaction provided that the jet is not Poynting flux dominated, 
i.e. the ratio of energy flux carried by magnetic fields and the kinetic en- 
ergy of matter {a) is not much larger than 1 . The mechanism described in 
this paper can also be applied to reverse shock emission. However, the ex- 
ternal reserve-shock has difficulties explaining the observed flat radio/mm 
lightcurves. 

* X-rays observed by the Swift satellite must have passed through the ex- 
ternal shock region as long as they are produced at the same, or smaller, 
radius as the external shock. The very high variabiHty of the x-ray data - 
on timescale of order minutes — suggests that the x-ray radiation is likely 
produced at a distance much smaller than the extemal shock radius. 
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Figure 1. Shown here is a schematic sketch of a jet produced in TDE. The 
energy is converted into x-rays at some distance from the black-hole by 
a as yet unknown mechanism. A part of this x-ray radiation is scattered 
by electrons in the forward-external shock which results in their efficient 
cooUng at early times. 



Finally, we calculate electron distribution function when IC cooling 
due to an external x-ray radiation field occurs in the K-N regime. 



3.1 Importance of IC cooling and its effect on electron 
distribution 

Let us consider a x-ray radiation source with luminosity Lx{tobs) 
(isotropic equivalent) in frequency band 1^1-1^2 which is produced 
by a source of size r^- Its specific luminosity scales as L^ oc v~^; 
all unprimed variables are quantities measured in the rest frame of 
the host galaxy of the source. Now consider electrons in external 
forward shock at a radius R(tobs), and the Lorentz factor of the 
shocked plasma (wrt CNM rest frame) to be T{tobs)- A photon of 
frequency u moving at an angle 9 wrt the radial direction is Doppler 
shifted to frequency i/' in shock comoving frame: 



u' ^,yr{i~pcose)^^[i + 6^r^], 



2r 



where 



.-2\l/2 



(1) 



(2) 



and we used Taylor expansion for cos 9, and F ^ 1 to arrive at the 
approximate expression for v' . From hereon all primed variables 
denote shock comoving frame. 

The specific intensity in the host galaxy rest frame is 



/.(e) 



L, 

4-K^rl 



(3) 



for 9<rx/R = 9x, and is zero otherwise. The specific intensity 
in the comoving frame of the shocked fluid is obtained by an ap- 
propriate Lorentz transformation, e.g. Rybicki & Lightman (1985), 
and is given by 

/,:, =/.(z.'H^~/./(8r^), (4) 

when 9 <C F^^. Therefore, the specific flux in the comoving frame 
is 



K' ~ {■^s'x)K' 



87rFi?2 ' 



(5) 



where 9'^ ~ 2r9x <S 1 is the angular size of the x-ray source as 
seen in the shock comoving frame. And the frequency integrated 
flux is 



J X 



(6) 



m-KV^R'^' 

The factor F^ in the denominator is easy to understand as it is due 
to Lorentz transformation of energy density. 

The equation for IC cooling of an electron of Lorentz factor 



7elS 



d(meC^7e) 



dt' 



4 „, ,2 



(7) 



where ar is Thomson scattering cross section. At the moment we 
are ignoring K-N effect which we will take up shortly. Thus, IC 
cooling time is given by 



t- 



(6 X 10''s)i??7L-^47(rV7e), 



(8) 



(^Tfxj'e 

where we made equation of equation ^ for /^, and have adopted 
the widely used notation X„ = X/10". The IC cooling calcu- 
lation assumed that the external shock region is optically thin for 
x-ray photons. This condition is easily satisfied since the optical 
depth to Thomson scattering is ~ arrieR ^ 10~^ UeRn ^ 1 
as long as Ue, the particle density in the unshocked CNM, is less 
than lO^cm"^ at a distance of lO^^cm from the blackhole which is 
indeed the case as shown by radio observations. Also, absorption 
of x-ray photons by the inverse-Synchrotron process is negligibly 
small. 

The dynamical time in the comoving frame is t'^^ ~ i?/(cF). 
Therefore, 



t' 



dy 



2VR17 

■y'eLx,4r ' 



(9) 



The Lorentz factor of the external shock a week after the tidal dis- 
ruption is about 2 (see §3.2 for shock dynamics), Rir ~ 1, and 
Lx ^ 5xl0'*'^erg/s. Thus we see from the above equation that the 
IC cooling time for electrons with 7^ > 30 is shorter than the dy- 
namical time. Radio and mm photons observed from Sw J1644-I-57 
are produced via synchrotron process by electrons with 7^ > 50, 
and so IC cooling cannot be ignored in the calculation of external 
shock radiation from this event. It should be pointed out that the 
synchrotron cooling is less important than the IC cooling at least 
for those electrons that are radiating well below the infrared banqj. 
For electrons with 7^ < 10 the IC scatterings take place in the 
Thomson regime for Sw J1644-I-57. The reason is that the energy 
of a 10 keV photon — Swift/XRT energy band is 0.3-10 keV — 
in the electron comoving frame is (107e/r) keV which is less than 
rrieC^ as long as 7^ < 10^. Higher energy electrons scatter 10 keV 
photons with a reduced cross-sectior|f|, and the scattered photon en- 
ergy is less than it is in the Thomson regime, i.e. < u'j'^. These 
two effects reduce the IC cooling of electrons with 7^ -J 10^ . A sim- 
ple, and fairly accurate, way to include these effects in IC cooling 
calculations is by only considering energy density in radiation field 
up to a photon energy of ej.„ : 



^kn 



2 / / 

rrieC /7e 



^kn ^knl ^1 



(10) 



^ The ratio of IC and synchrotron cooling times, in Thomson scattering 
regime, is equal to the ratio of energy density is magnetic field and radiation 
field. The x-ray luminosity from Sw J 1644+57 was very high for the first 
several weeks, and thus it turns out that the energy density in radiation is 
much larger than magnetic energy density as long as eg for the shocked 
fluid is smaller than 1 . 

^ Klein-Nishina cross-section is smaller than the Thomson cross section by 
a factor of photon energy in electron rest frame divided by rrieC^. 
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where h is Planck's constant. The spectral index in the x-ray band 
for Swift J1644+57 is /3 ~ 0.7, and in that case the reduction to 
the IC cooling due to K-N effect scales as 7^^ ' for large 7^ (see 
§3.1.1). This is a weak effect but has important implications for 
electron distribution function which is discussed below. 



3. 1. 1 Electron distribution function in external shock 

Electron distribution function, dne/d'ye, in a shocked plasma is de- 
termined by the dual effect of acceleration at the shock front and IC 
and synchrotron coolings while electrons travel down-stream. The 
distribution function is obtained from the following equation: 



d dUe d 
dt d'ye d'ye 

where 

4(tt 



"fe 



drie 

d'ye 



= S(7.), 



7.! = 



SnieC^ 






7e/3e, 



(11) 



(12) 



is the rate of change of 7e due to IC and synchrotron losses, /(< 
Vkn)lc is energy density in radiation below the K-N frequency i^fc„ 
(eq.[T0l(, and 



Side) CCye"^ for 7e ^ 7i 



(13) 



is the rate at which electrons with LF 7^ are injected into the sys- 
tem; 7i is the minimum Lorentz factor for shock accelerated elec- 
trons. All the variables considered in this subsection are in shock 
comoving frame. However, we will NOT use prime (') on variables 
in this sub-section (and only in this sub-section) to denote shock 
comoving frame as this is the only frame being considered here. 
When results from this sub-section are used elsewhere in the paper 
we will revert back to the proper notation for prime and un-prime 
frames. 

Let us define radiative cooling time 



tc(7.) = ^, 

and cooling Lorentz factor, 7c, which is such that 

ic(7c) = I'dy, 



(14) 



(15) 



where tdy is the dynamical time. 

We will calculate the distribution function in two cases. The 
slow cooling case where 7^ > 7i, and the fast cooling regime when 
7c < 7i. Moreover, we will consider here only IC losses in K-N 
regime as that is relevant for the "afterglow" radiation of the TDE 
Swift 1644-1-57. We note that the theoretical lightcurves and spectra 
shown for this event in the next section are obtained by numerical 
solutions of relevant equations that include synchrotron loss and no 
assumption regarding K-N regime. 

Since vu-a oc 7^^ (eq.fTO)) and /(< Ukn) oc 7e 

oc 7c^ . We see from equation dllt that for 7e < jt 

1 



-(1-/3) 



there- 



fore 7e 

dUe 



d'ye Je 

whereas for 7, 
dn, 



>7i 



I,:d^Si'y) 



(16) 



(17) 



d-Je ie 

Therefore, the solution of equation ( 111) in the slow cooling case is 



dne 

d-ye 



-P-/3 



7i ^ 7e ^ 7c 
7e >7c 



(18) 



and in the fast cooling case 



dUe 

d'ye 



-1-/3 
-P-/3 



7c ^ 7e < 7i 

7c > 7i 



(19) 



The synchrotron spectrum for particle distribution dUe /d'ye oc 



7c ' is fv 



-{9-l)/2 



for 1/ > ramive, Vi), and it is v 



1/3 



for 



V < nim{i'e,Vi); these spectra are modified when synchrotron- 
absorption becomes important (the criterion for this is discussed in 
the next sub-sectionjj. 

It should be clear from these equations that the spectral in- 
dex above the peak depends on p, and the spectral index /3 of the 
radiation field that is responsible for IC cooling of electrons; the 
spectral index for i^ > niax(i/i, Ue) is (p — 1 -I- /3)/2 which is 
harder than v~^''^ . As an application to Swift J1644-I-57 we note 
that the spectral index in the radio band above the peak of the spec- 
trum is reported to be ~ 0.75 for fobs ^ lO^d (Berger et al. 2012). 
This does not mean that p = 2.5 as suggested by eg. Berger et al. 
If IC cooling in K-N regime is important for determining particle 
distribution, which equation (|9) for t'^^/t'^y makes clear is the case, 
then p = 2.3, since according to the late time X-ray data (3 = 0.4 
(Saxtonetal. 2012). 



3.2 Shock dynamics and synchrotron lightcurves 

The dynamics of external forward shock can be determined, ap- 
proximately, from the following energy conservation equation 

471 



(3-.) 



i?^-Ri7n/mpC^r(r -l)= E, 



(20) 



where E is the energy (isotropic equivalent) in the blast wave, the 
density of the CNM is taken to be 



%{R) = -nfR^Y . 



(21) 



and Rn — _R/10^^cm. For a relativistic shock, i.e. F 2> 1, the 
dynamics is described by 



(3-s)£ 



'i-irrnpC^ri f R^ R^j 



1/2 



(22) 



and the deceleration radius - where the relativistic ejecta has im- 
parted roughly half its energy to the surrounding medium - is 

(3-s)£51 



R. 



'inmpC'^rifT'^ ' 



(23) 



where Fo is the initial Lorentz factor of the relativistic outflow. For 
a non-relativistic shock the solution of equation ( 120) is 



(3-s)£51 



271771,^71/ _Rj7 ' 



1/2 



(24) 



The Lorentz factor for the external shock of TDE Swift 1644-1-57 
was of order a few at early times (less than a week), and the shock 
speed decreased to Newtonian regime at late times. Hence the an- 
alytical calculations we present in this section considers both of 
these regimes. The numerical results described in the next section 
are more accurate and do not rely on making any assumption re- 
garding F. 

A photon released at radius R arrives at the observer at time 



^ A very detailed description of synchrotron spectrum when IC cooling in 
K-N regime is important can be found in Nakar et al. (2009) 
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tob 



(1 + 2) 



dR 



1 



(25) 



For relativistic blast waves tots oc R^'", and for non-relativistic 
case tabs oc w''~"'''^. Moreover, the shock speed as viewed by an 
observer is 



roct„ 



and 



Roct 



l/(4-s) 
obs ' 



for the relativistic case and 

3-3 



w oc t 



and 



Roct 



2/(5-3) 
obs ' 



(26) 



(27) 



for Newtonian dynamics. 

The deceleration time in observer frame, for relativistic out- 
flow, is given by 



fd « (3 X 10*s) 



(3 - s)E5 



1.8n 



f,2 



(i + 2)ro7 



(28) 



Thus, for E = lO^^erg, n/ = lO^cm"^, and To = 10, the de- 
celeration time is of order 5 hours, and the Lorentz factor of the 
shocked CNM 10 days after the start of TDE is ~ 2. 

3.2.1 External shock radiation and lightcurves 

We provide analytical calculation for synchrotron radiation and 
lightcurves from external forward shock in this sub-section; more 
accurate numerical results are presented in the next section. 

The magnetic field strength and the electron minimum Lorentz 
factor in a shock heated plasma are estimated by assuming that a 
fraction of the energy of shocked plasma (eg) goes into the genera- 
tion of magnetic fields, and a fraction (ee) is taken up by electrono 
Under this assumption we find that 

B' = [327reflmpC^n(7?)r(r - 1)] ^^^ , (29) 

and the minimum Lorentz factor of shock accelerated electrons is 

, _ (p - 2)mp 



-^e(r-i) 



(30) 



{p - l)me 
Using results for shock dynamics from last sub-section we find 



B' oc 



t„ 



t„ 



r/3> 1 



r/3< 1 



(31) 



and the synchrotron frequency corresponding to electron Lorentz 
factor 7,' is 



"7fr 



t 



'3/2 



2-KmeC 



r/3>i 



r/3< 1 



(32) 



We next calculate synchrotron frequency for those electrons 
that cool on a dynamical timescale (uc) due to IC scattering of x- 
ray photons. Let us rewrite equation (|7) to explicitly consider x-ray 
flux below the K-N frequency 



dmeC^7e 
dt'„ 






(33) 



° This crade approximations is made due to our inability to calculate mag- 
netic fields and electron acceleration in collisionless shocks ab initio. And 
the dimensionless parameters eg and te hide our ignorance regarding these 
processes. 



where u'^^ is given by equation l IlOl ). The LF of electrons that cool 
on a dynamical time is obtained from the above equation 



SrUeC 



For a powerlaw spectrum for x-ray radiation 

Therefore, the electron "cooling" Lorentz factor is given by 



7c oc 



2pl+/3 



R-'V 



1//3 



(34) 



(35) 



(36) 



fdy'^^V'oi'') 

Using equations for shock dynamics from last sub-section we find 

ni//3 



7c oc < 





^[2-(2 + ,3)(3-5 


)l/[2(4- 


= )1 








ix(t„6.) 






■,(»-l)/(5- = )- 


1//3 




ix(t„6.) 





for r/3 > 1 



for r/3 <1 



(37) 



The synchrotron frequency in observer frame corresponding to 7^ 
is 



?'7c'r 



27r(l + z)Tn,^c^ 



(38) 



Or 



I^c oc < 



-3/2 
'■obs 

,-3/(5 


f*ir 


-2)/(i-s)' 


2/(3 


L 


x(to6s) 
■tCs-l)/(5 


- = )■ 


''obs 


ix(tob. 


) 



2//9 



for r/3 > 1 



for r/3 < 1 



(39) 



For s — 1.5, /3 — 0.7 and Lx oc t^^^^ wc imu i^c 'j^ i-obs 



-5/3 



we find fc oc t^j^ for 
relativistic case, and i/c oc t^ji, for Newtonian dynamics. For these 
parameters for Sw J1644-I-57, Vc is about 1 GHz at 10 days, and it 
rises to 100 GHz at ~ 70 days. 

The observed flux at a frequency v > niax(i/a , i^i , J'c ) is given 
by 



where Ua is synchrotron absorption frequency. 



/p^ 



q'^B'TNe 1 + z 
rrieC^ 47rd|^ ' 



(40) 



(41) 



TVe OC R^~^ is the total number of electrons with Lorentz factor 
> min(7j', 7^), and d^ is the luminosity distance to the source; we 
made use of equations J18t and M9\ in deriving equation ( 140) . 

Making use of equations ( |26t , Ml\ , ( 1311 ), ( |32| ( and J39b in the 
above expression for /^ we find the observed lightcurve behavior: 



fu{tobs) oc 



-(p+/9-l)/2 



i-'x [Job 



-[2+3/3+3(p-l)]/4 



19-63-p(15-4B)-3ff 



[2(5-s)l 



^ obs 



r/3>i 



r/3<i 



(42) 



-5/3 



Forp = 2.0, P = 0.7, s = 1.5 and Lx oc f "j/ , fu ^ ^ v^^^ 

for relativistic dynamics, and /^ oc //""^^fo^j for Newtonian 
shock. The dynamics of Swift J1644-I-57 for lOd < tobs < lO^d lies 
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between these two extremes, and the lightcurve is flat for as long a 
duration as the observing stays above Vc- Of course, at sufficiently 
high frequencies — which of order 10^'^Hz for this TDE at ^ lOd 
— electron cooling by IC scattering becomes unimportant due to 
K-N suppression, and synchrotron cooling time becomes smaller 
than the dynamical time. In this case the lightcurve decline has the 
well known form of /^ oc !/^''"f~j^\j ''" '' for relativistic shocks, 
e.g. Kumar (2000). 

The rapid increase of Uc with time is responsible for the nearly 
flat lightcurve observed for Swift J1644+57 in frequency bands 
i^ > max(!/a, Vi, Vc)- When fcitobs) moves above the observing 
band then the lightcurve resumes its usual, ~ f"^,^, decline which 
occurs at tots ^ lO'^days for v > lO^GHz as shown in numerically 
calculated lightcurves presented in the next section. 

Another interesting case to consider is when Vc < i^ < i^i and 
the observing band is above the synchrotron absorption frequency. 
The lightcurve in this case is determined from 



J- /, \ /• -/3/2 8/2 

Making use of equations l l41b and l |39l > we find 



jv(toba] OC 



' -(2+3/3)/4 
-/3/2 ^obs 



Lxitobs) I [4-2s-3/3]/[2(5-3)] 
''obs 



r/3>i 



r/?<i 



(43) 



(44) 



For , 



0.7, s = 1.5 and Lx oc t 



-5/3 



''obs '■^^ lightcurve rises as 
j^o.35 ^0.64 |.^i^5j fpj. j-elativistic (Newtonian) blast waves. 

Next, we calculate the synchrotron self-absorption frequency 
(Va) when IC cooling is important. The easiest way to determine 
fa is by using Kirchhoff's law or the fact that the emergent flux at 
frequency Va is approximately equal to the black-body flux corre- 
sponding to electron temperature given by 



kTl, = max[7l,min(7i,7^)]meC^ 



(45) 



where 7^ is electron Lorentz factor corresponding to synchrotron 
frequency u'a, i.e. j'a oc {v'a/ B'Y'^ . The equation for v'a, using 
Kirchhoff's law, is 



2kTy^ aTB'm,&N^{> kT^/m^cf) 



Anq 



(46) 



where the left side of this equation is the Planck function in the 
Rayleigh- Jeans limit, and Ne{> kTe/nieC?) is the column density 
of electrons with LF larger than max[7o, min(7i, 7^)]. 

One case of particular interest for application to Sw J1644-I-57 
is where Va is larger than inax(!/i, i^c). We can rewrite equation 
l|46t for this case as 



2mt 



/5/2 



/ 27r7ne 



V qB 






(47) 



Making use of equations | |26I ). l l27b . Oil ) and l l42t we find 



i^a cc < 



2(4-a)(p + l3+i) 

bs 

2 



30-6s-p(15-45)-3/3 

(5-s)(p + ;3 + 4) 



P + fi+i 



r/3>i 



r/3<i 



(48) 



5/3 



For s = 1.5, p ^ 2, 13 = 0.7, and L^ oc t~^J", Va oc t^^^ 
i^obs) for relativistic (Newtonian) shock dynamics. We note that 
due to the rapid rise of Vc (see eq. 1391 ) the above equation for Va 
ceases to be applicable after a few lOs of days for Sw JI644-I-57 



when the ordering of synchrotron characteristic frequency changes 
to Vi < Ua < i^c- The synchrotron absorption frequency for this 
case can be shown to be 



Va oc 



1 
P+4 


10-s 
2(4-s) 


, 3(p- 
+ T 


ly 


obs 








[(5+4 


) + (15 


-4s)(p- 


1)1 




(5~3)(p+4) 





r/3»i 
r/3<i 



(49) 



For s = 1.5 and p = 2,VaOC t~°/^^ (^Jb°/^) for relativistic (New- 
tonian) shock dynamics. 

The rise of lightcurve for about 10^ days for u < lOGHz is due 
to the fact that these observing bands lie below the synchrotron 
absorption frequency during this period. The observed flux below 
i^a decreases rapidly with decreasing i^, and hence even a slow de- 
crease of Va with time leads to a rising lightcurve. 

We note that the luminosity in high energy photons produced 
by the IC scattering of x-rays by electrons in the external shock is 
of order L^Y, where Y ~ aTn{R)R'y'^-y[ ~ 10~^-Ri7n/_2 is the 
Compton parameter, and the energy of these scattered photons is 
~ lOMeV. 



4 NUMERICAL RESULTS AND APPLICATION TO 
SWIFT J1644+57 

Equations l l20t & l l25t are solved numerically to determine shock 
dynamics as a function of observer time, and equations JIl) . il2t . 
l |29l l, l |30l l & l |46t are used for calculating spectra and lightcurves 
that take into account electron cooling due to IC scattering of an 
"external" radiation field with flux and spectrum taken from x-ray 
observations, and also cooling due to synchrotron radiation. No as- 
sumption regarding whether IC scatterings are in K-N or Thomson 
regime are made for calculations presented in this section. 

The energy in the blast wave is injected only at the beginning 
of numerical calculation when the blast wave is at a radius R — 
lO^^cm; results are insensitive to the precise value of the initial 
radius as long as it is taken to be much smaller than the deceleration 
radius. In Fig. 2 we show R, F, B' ,^i & 7c, for parameters relevant 
forSwJ1644-H57. 

Theoretically computed lightcurves for Sw J1644-I-57 in fre- 
quency bands 230 GHz & 87 GHz, together with observed data 
from Berger et al. (2012), are shown in Fig. 3. The theoretical 
lightcurves were calculated from a decelerating blastwave which 
had a fixed amount of total energy (i.e. no late time injection of en- 
ergy) and yet the observed flux in these frequency bands is nearly 
flat for ~ lO^days (consistent with the observed data) in spite of 
the fact that these bands clearly lie above the synchrotron injection 
(Vi) and self-absorption (i/a) frequencies (shown in fig. 4). There is 
a noticeable dip in lightcurves between ~ 10 & 20 days. The reason 
for this dip is that the observed x-ray lightcurve during this period 
is nearly flat (except, of course, for the ever present fluctuations). 
A flat x-ray lightcurve means that Vc does not increase during this 
time interval (see fig. 4). However, Vi and the peak flux (which are 
independent of Lx) continue their monotonic decline, and that is 
the reason for the flux decrease during this period. A more realis- 
tic treatment of the x-ray flux in our calculation, that includes the 
observed fluctuations, should provide a better fit to the observed 
data. 

The nearly flat lightcurves for ~ 10^ days for 87 and 230 GHz 
is due to the fact that these bands lie above the synchrotron cooling 
frequency which rises very rapidly with time (VcZ tabs) due to the 
rapidly declining x-ray flux; we obtained this behavior analytically 
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Figure 2. A few of the parameters of the external forward show are 
shown in this figure as a function of observer time. Shock front radius 
(top left), Lorentz factor of the shocked fluid (top right), magnetic field 
in the comoving frame of shocked fluid (bottom left), and bottom right 
panel shows 7^ (solid line) & 7^ (dotted line). The energy in the blast- 
wave was taken to be B = 1.5 X 10^'^erg (isotropic), the CNM den- 
sity n{R) = (400cm-3 X R~^, eg = 4.5 X IQ-^, ee = 0.1, and 
p = 2.3 for these calculations. The x-ray luminosity was taken to be 
L^ = 4 X 10'"'erg/s(toi,s/20days)-5/3 for tobs > 20d and constant 
between 10 and 20 days; x-ray photons passing through the external shock 
cooled electrons due to IC scatterings and that was the dominant cooling 
mechanism for electrons radiating at radio/mm frequencies for ~ 10^ days. 
We note that the uncertainty in 7c is small as it depends primarily on Lx 
which is well determined from observations. However, ■ji might be much 
larger than shown in this figure, particularly for p Ki 2 favored by our 
modeling of the observed data, if the minimum Lorentz factor of electrons 
accelerated in shocks is larger than the value we have adopted in this paper, 
viz. 7i = Eeip - 2)(r - l)nip/[me(p - 1)]. 



in the last section, and the numerical result is presented in fig. 4. 
The lightcurve remains nearly flat as long as the observing band 
is above Uc, and then it starts to decline as ~ i~. ' . The duration 
of the flat phase is dictated by the temporal behavior of x-ray flux 

— 5/3 

which we took for these calculations to decline as t^^^ , and the 
spectrum was taken to be a powerlaw function; /^ oc i/""'^ for 
O.lkeV < i^ ^ lOkeV as per Swift/XRT observation^ for tots < 
40 days, eg. Saxton et al. (2012). The observed x-ray lightcurve 
in reality is much more complicated, and shows rapid fluctuations. 
However, we don't expect Vc, and radio lightcurves, to follow these 
wild fluctuations because the signal we receive from the external 
shock at any given fixed observer time is in fact an integral over 
equal arrival time 3-D hypersurface that has a large radial width so 
that the effect of x-ray radiation on electron distribution function is 



^ Swift/XRT observes photons with energy between 0.3 and 10 keV. How- 
ever, it would be a remai'kable coincidence if the spectmm for Sw 71644-^57 
starts deviating from the observed powerlaw behavior exactly at the low or 
the high energy threshold for XRT. The infrared k-band flux at 10 days was 
0.03 mJy (Levan et al. 201 1) and at the same time the mean 1 keV flux was 
0.02 mJy (Zauderer et al. 2012), and that suggests that the XRT spectrum 
became much harder, perhaps u^'^, at ~ 0.1 keV. The flux above 10 keV, 
however, is rather poorly constrained by observations. 
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Figure 3. Lightcurves in frequency bands 230 GHz (solid line) & 87 GHz 
(dotted curve) calculated for the model described in §3, together with ob- 
served data for Sw J1644-H57 from Berger et al. (2012) are shown in this 
figure. The pai'ameters used for these calculations are same as that in Fig. 2 
(see fig. 2 caption for details). Lightcuree before 10 days is not shown here 
because the wild fluctuation in x-ray flux at earlier times (by more than a 
factor 10) introduce huge error in our theoretical calculations which is based 
on a powerlaw decline model for x-ray luminosity for tots > 20days, and 
a constant luminosity for lOd ^ toj^ ^ 20d. 



averaged over a time duration Stots of order tots- The x-ray flux 
averaged over this time scale is not a smooth powerlaw function, 
but consists of periods of faster and slower declines. The duration 
of time when lightcurve in an observed frequency band is roughly 
flat depends on these fluctuations, and for that reason the duration 
of the flat phase shown in figure 3 should be considered to have an 
uncertainty of order 20%. 

The lightcurve decline for frequency i' which lies below the 
synchrotron cooling frequency (vc) is easy to obtain using the 
scalings provided in the §3.2.1. For sub-relativistic shocks, B' oc 



t 



-3/(5-s) 



(eq.[3Tjl, and ji oc i^^ 



(15-4s)/{5-s) 



(eq. |32t , and so the 



observed flux for Uc > u > ma.x{i'a, Vi) decreases with time as 



-[{15-4a)(p-l)-(6-4s)]/2(5-s). 



for J1644-I-57, Uc ~ 10 GHz at 
tabs ~ lO^d (Fig. 4). For s = 2 and p = 2.3, the lightcurve de- 
cline is expected to be t"^^' which is consistent with the observed 
decline (see. the 87 GHz lightcurve in Fig. 3 on this page, and also 
8.4 — 43 GHz lightcurves in Fig. 3 of Zauderer et al. 2012). 

The reason that we don't show theoretical lightcurves in fig. 
3 for tabs < 10 days is that the x-ray lightcurve before 10 days is 
dominated by very large flares — the flux decreased by a factor '^ 
10^ between 4 and 6 days — which cannot be captured by a simple 
powerlaw model used in this paper. A calculation of IC cooling of 
electrons, and the emergent flux, before 10 days must include these 
wild fluctuations for the result to be meaningful. 

The rise in lightcurves at lower frequencies (10 GHz or less) 
and for lOd < toba ^ lO^d is due to Va lying above the band (see 
Fig. 4). This behavior can be understood using analytical scalings 
presented in the last section. The above remark regarding fluctu- 
ations in Lx{tobs) applies to Va as well, although self-absorption 
frequency is much less sensitive to L^ than Vc is. 

The solutions we are finding has min{vi, i^c) < Va, and that 
might naively suggest that the spectrum below Va should be v'^'^ 
in conflict with the observed spectrum of ~ v^'^ for tabs ^ lO^d. 
The fact that the observed spectrum is shallower than v^ is already 
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Figure 4. The left panel shows the three synchrotron characteristic frequen- 
cies, obtained by numerical calculations, for the same set of parameters as 
fig. 2; see fig. 2 caption for details; solid line is i/^, dotted line is Vi and the 



dashed curve is v^. Note the sharp rise in i^c, as ~ t'^ 



, between 20 and 100 



days, which is consistent with our analytical estimates (see §3). We note that 
the uncertainty in the calculation of v^ is small. However, the actual value 
of Ui might be larger than in this figure by 1-2 orders of magnitude for 
p Si 2 favored by our modeling of the observed data (see caption for Fig. 2 
for explanation); good fits to observed lightcurves (like those shown in fig. 
3) can be obtained even when Ui is much larger than shown in this figure by 
appropriately rescaling microphysics parameters and changing E and n by 
a factor ^ 2. The right side panel shows the spectrum calculated at 20 days. 
The spectrum above Vi and u^ is v~^-^ even though we took p = 2.3. 
This is due to the fact that when cooling is dominated by IC scatterings of 
x-ray photons in the K-N regime, the spectrum is ;y^(P^i+'9)/2 instead of 
y— p/2. ([jg spectral index for x-ray radiation was j3 = 0.4 for tabs ^ 50d 
(Saxton et al. 2012). The spectrum steepens to ~ v~^-^^ for v > lO^^Hz 
when the cooling is dominated by synchrotron loss. 



below the peak. The reason is that when electrons are cooled by IC 
scatterings in K-N regime the resulting spectrum is y~(P~^+°'i'^ 
(see §3.1) and not v^'^''^ . X-ray data show that /3 varies with time 
(Burrows et al. 2011), and that /3 was about 0.7 for toba ^ 50d and 
~ 0.4 for t > 50d (Saxton et al. 2012). Therefore, we expect the 
radio spectrum above the peak to be ^ i/""*^ at early times, and 
i^""'^ at late times, for p = 2. The late time radio data are con- 
sistent with this expectation, and for fobs ^ 50d spectra above the 
peak are poorly sampled. 



4.1 High energy radiation 

IC scattering of an x-ray photon of frequency V:^ by an electron 
of Lorentz factor 7^ in the external shock produces a high energy 
photon of energy ~ 2'^'I^Vx. Considering that 7^ ~ 30 (Fig. 2), the 
IC photon energy is ~ 10 MeV, and the expected IC luminosity is 
~ ttLi:7- min(7c,7i); where tt ~ aTn{R)R ~ 10~^i?i7W/,2 
is Thomson optical depth. At tabs = lOd, Lx ^ 3 x 10'*® erg 
s~^> 7c ^ 'iRn & ji ^ 30 (fig. 2), and therefore we expect 
IC luminosity to be ~ 5 x W^'^Rirnf^2 erg s~^. Another way 
to estimate the luminosity of high energy IC photons follows di- 
rectly from the fact that IC scatterings of x-ray photons is effi- 
cient at cooling electrons at early times, i.e. IC cooling time is 
less than or of order the dynamical time for tabs < 20d. Hence, 
the IC luminosity should be roughly equal to the energy lumi- 
nosity carried by electrons accelerated by the external shock, i.e 



L^c ^ 16TTeeR^mpC-'n{R)r^ 



-in/,2-R?7erg s ^. 



There is no data between 1 & 100 MeV for Sw J 1644-1-57 to our 
knowledge, however, Fermi/LAT provided an upper limit of about 
W'^erg/s above 100 MeV at lOd (Burrows et al. 2011). This sug- 
gests that 71/, the CNM density at i? = lO^^cm, cannot be larger 
than about lO^cm"^. 



telling us that something is missing. And most likely what that is, 
is that the jet has some angular structure, and when integration 
over equal arrival time hyper-surface is carried out to obtain the 
observed spectrum, it turns out to be softer by u"'^ — f^ than the 
comoving, instantaneous, spectrurn^'^L 

The infrared flux at 10 days (0.03 mjy, Levan et al. 2011) 
is consistent with the theoretically calculated value for parameters 
used for fitting 230 GHz lightcurve shown in Fig. 3. The infrared 
producing electrons have 7^ ~ 10* and their IC cooling by scat- 
tering x-rays photons is highly suppressed due to K-N effect. How- 
ever, Vc due to synchrotron cooling is ~ 10^'^Hz, so the spectrum 
in this band should be ~ 1/^^. Levan et al. (2011) report that the 
infrared lightcurve declined with time as ~ t~^^ between 5 and 10 
days. This is much too steep for external shock and we suspect that 
this might have something to do with the shaip drop in the x-ray 
flux between 4 and 7 days by a factor ~ 50. 

A theoretically calculated spectrum for external shock syn- 
chrotron radiation at 52 days is shown in Fig. 4. The spectrum 



above the peak is quite hard. 



even though Ui and Uc lie 



^^ The equal arrival time hyper-surface is a 3-D section of the 4-D space- 
time with the property that radiation from every point on this hyper-surface 
arrives at the observer at the same time. The radius (R), 7^, B' etc. for dif- 
ferent points on the hypersurface are in general different for a non-uniform 
jet, i.e. a jet where energy density, Lorentz factor etc. vary with angle wrt 
jet axis and distance from the jet surface. It can be shown that an integra- 
tion over equal arrival time hyper-suiface invariably softens the observed 
spectrum. 



4.2 Parameters of the TDE and the external shock 

An accurate determination of jet energy and other parameters for 
Sw J 1644-1-57 should take into consideration the effect of fluctuat- 
ing Lx on electron distribution, a more detailed treatment of the 
dynamics than described by equation ( |20t , and the 3D structure of 
shocked plasma and magnetic fields. But that is beyond the scope 
of this worH^^l We show in Fig. 5 an approximate determination 
of E (isotropic energy in the external shock), the CNM density, 
and e_B (energy fraction in magnetic fields) allowed by the multi- 
wavelength data for Sw J 1644-1-57; see the caption for Fig. 5 for a 
description of the various observational constraints that were used 
for the determination of these parameters. 

Fermi/LAT upper limit restricts the CNM density at lO^'^cm 
to be smaller than ~ 5x10"^ cm^'' (see ii4.lt . And according to 
the radio/mm/infrared data the isotropic energy in the blastwave 
is '^ 2 X 10^"' erg which is consistent with estimates obtained using 
equipartition argument, e.g. Barniol Duran et al. (2013). We note 
that the uncertainty in these parameters and eb is a factor of a few. 
It should be pointed out that the afterglow data we have used for 
determining energy etc. consists of observations carried out 10 or 
more days after the TDE. The blastwave Lorentz factor at this time 
had dropped to < 2, and the angular size of the shocked CNM is 



The goal of this work is to show that cooling of electrons by x-ray 
photons passing through the external shock can explain the puzzling, flat, 
lightcurves observed in mm/cm bands for ~ lO^days. 
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Figure 5. Isotropic blast wave energy (E), CNM density at R = 10^^ 
cm, and eg allowed by the data for Sw J1644+57. The parameters shown 
in this figure should be considered a rough estimate due to the approxi- 
mate way we have handled the blastwave dynamics, and used a smooth 
decline for the x-ray flux passing through the external shock and cooling 
electrons instead of the highly fluctuating flux Swift/XRT observed. The 
observational constrains that were used in the determination of parame- 
ters are: (1) flux at 230 GHz at 10 days after the tidal disruption was 10 
mJy (Berger et al. 2012), (2) Ua = 15 GHz at 10 days, (3) the flux at the 
peak of the spectrum at 10"^ days is 50 mJy, (4) an upper limit of 10^® erg 
s~ for >10 MeV at 10 days which suggests that nf < 10 cm~ (rif is 
density at R = lO^^cm); a tolerance of ±30% was used for all these con- 
straints which reflects the degree of uncertainty of our calculations. We took 
Lx = 4x 10''®erg/s(to6s/20days)-5/3 for tabs > 20d and constant 
between 10 and 20 days, and e^ =0.1 for these calculations; for e,. = 0.5 
these plots look very similar, except that they are shifted to the left by a 
factor ~ 10. 



perhaps larger by a factor ~ 2 than the jet angle at the smaller ra- 
dius where x-rays were produced. Therefore, the beam corrected 
energy in the blastwave is perhaps a factor ~ 2-4 larger than the 
energy in x-rays which suggests a > 20% efficiency for converting 
jet luminosity to x-rays. 

We have also looked for the density stratification (s) and find 
that s > 1.5 provides a reasonable fit to the available data; s = 1.5 
density stratification is expected for an advection dominated accre- 
tion flow in the neighborhood of a blackhole, and s = 2 is when 
outflow determines the CNM density structure. 

Although we took e^ =0.1 for results shown in all of the 
figures, other values of e^ are allowed by the data as long as we 
takeS-e-\ 

We find that flux in the infrared K-band at 10 days for the 
solution space shown in Fig. 5 is between 0.02mJy and 0.04 mJy; 
the observed K-band flux at this time was 0.03 mJy (the error in 
flux calculations, due to various simplifications we have made, is a 
factor < 2). The observed K-band lightcurve for tabs 5; lOd decayed 
as ~ t~^^ (Levan et al. 201 1) which is too steep for external forward 
shock emission and we suspect that this might have something to 
do with large x-ray flares that occurred between 5 and 10 days (see 
§4.3 for further discussion). 

The 1 keV x-ray flux at 610 days after Swift/BAT trigger — 
for the entire solution space shown in Fig. 5 — is within a factor ~ 
2 of the value observed by the Chandra satellite which was 0.7nJy 
(Levan & Tanvir, 2012). 



4.3 A puzzling drop in Va between 5 & 10 days 

According to Berger et al. (2012) the synchrotron-absorption fre- 
quency (ua) dropped by a factor 10 between 5 and 10 days after the 
tidal disruption event Sw J1655-I-57. If these observations are cor- 
reco then this result essentially rules out any shock model, with 
or without continuous energy injection, for radio emission unless 
some mechanism external to the shock led to a drastic change to 
the electron energy and magnetic fields in shocked plasma between 
5 and 10 days. To understand this result we make use of equation 
(|46l> for u'„ : 

u'^'-fU'^'a) cc B'N4^ 7^ or /.„ (X r [B'^/^N^^, 7:)]°'* (50) 

where 'y'^iy'a) is the Lorentz factor of electrons that have charac- 
teristic synchrotron frequency v'a, and Ne{^ 7e) is the column 
density of electrons that radiate above v'^. The right hand side of 
the first equation is proportional to the synchrotron flux at i/'a in 
the comoving frame of the source when synchrotron absorption is 
ignored, and in the second equation we have considered the case 
where v'^ > min(i^^,i^c) (other cases are straightforward to con- 
sider). For any external shock model, with or without energy in- 
jection and arbitrary CNM density stratification, F, B' and N^ are 
smooth functions of time (see §3). Since F oc t~j^ ""'/y ~ ^1 (gq 
I26t , the most change we can expect for F between 5 and 10 days is 
a factor 2^/* ~ 1.3 for a s = medium; similarly, B' oc TR-''^^ 
can't change much during this time in the absence of an external 
agency at work. And A^e, which depends on the number of elec- 
trons in the CNM swept up by the shock, is typically an increasing 
function of time and hence its effect on u'^ is to increase ito We 
know empirically that Ne{^ l'e)B' oc f{ua)/V decreased by a 
factor ~ 3 between 5 and 10 days for Sw J1644-I-57 (see fig. 2 of 
Berger et al. 2012), and that by itself can only account for a de- 
crease of Pa by a factor < 2. We see, therefore, that unless there 
is a dramatic change in B' and/or "fi(p'a) on a short time scale of 
Stobs/tobs = 1 — which is not possible in external shocks without 
a powerful radiation front passing though it or some such mecha- 
nism — Va cannot drop by a factor ~ 10 between 5 and 10 days 
that is observed for Sw 1644-1-57; any addition of energy to a decel- 
erating shock only makes the problem worse. 

It is probably no coincidence that the x-ray luminosity of Sw 
1644-1-57 fell by a factor ^ 50 during the same time period that 
Va dropped by a factor 10. When a highly luminous x-ray front 
was passing through the external shock (tots ^ 6d), electrons in the 
region were cooled by IC scatterings of x-rays to a Lorentz factor 
of order unity in less than one dynamical time (see eq. |9) . The sharp 
drop in x-ray luminosity during 6 to 8 days by a factor 50 reduced 
the IC cooling of electrons and raised 7^ by the same factor. As we 
can see from equation ( I50t a rising electron temperature helps to 
lower Va , however, the effect on 7' (I'a ) is much smaller. Ultimately, 
a sharp drop in B' is needed to account for the observed decrease of 
Va as can be seen from the second part of equation ( 150) . The rapid 



t^ We would not be surprised if it were to turn out that the drop in Va has 
been overestimated by observers by a factor 2-3. The spectrum reported by 
Berger et al. (2012) at t^bs = lOd has no frequency coverage between 15 
and 80 GHz. It is possible that the actual peak was at ^^ 50GHz, in which 
case Ua at lOd has been underestimated by a factor ~ 3. 
^^ If a shock were to run into a dense cloud, that would cause Ua to in- 
crease. And when a shock suddenly encounters a zero density medium, 
the shocked plasma undergoes adiabatic expansion and then 7^, B', and 
Afg(JS u'a) decrease with time. It is easy to show that the net effect of adia- 
batic expansion is a slight decrease of Va in time interval tabs — "^ioba- 
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fall off of the infrared lightcurve mentioned in §2 might also be 
due to this decline of magnetic field. We suspect that a sharp drop 
in Lx leads to a smaller B' . One possible way this might occur is 
that while L^ is large, or jet power is high, the interface separating 
the shocked jet and the shocked CNM is Rayleigh-Taylor unstable 
and the turbulence associated with this instability produces strong 
magnetic fields. Currents driven in the shocked CNM by high L^ 
might be another way to generate close to equipartition magnetic 
field in the region where radio/mm emission is produced. However, 
a connection between L^ and magnetic field generation is just a 
speculation at this poinO. 



5 CONCLUDING REMARKS 

An interesting puzzle is posed by the radio/mm data for the tidal 
disruption event Swift J1644+57; lightcurves in frequency bands 
that lie above the synchrotron characteristic peak frequencies, i.e. 
V > IQ GHz, are roughly flat for lO^days rather than decline with 
time as ~ t~^^ as one would expect for radiation produced in the ex- 
ternal shock. Berger et al. (2012) suggested that the flat lightcurves 
are due to continued addition of energy to the external shock for 
10^ days. This requires about 10 times more energy to be added at 
late times than is seen in x-rays, and that is worrisome. We have 
described in this paper a simple model for flat lightcurves that does 
not require adding any extra energy to the decelerating shock wave. 
The flat lightcurves, instead, are a result of efficient cooling of elec- 
trons in the external shock by IC scatterings of x-ray photons that 
pass through the external shock region on their way to us (Fig. 1). 
The IC cooling becomes weaker with time as the x-ray flux declines 
and that is the origin of radio/mm lightcurve flatness. This model 
can account for the radio, milli-meter, infrared and > lO^MeV data 
(temporal decay as well as spectra). 

The angular size of the x-ray beam should not be much smaller 
than the size of the external shock region during the initial ~ 
lO'^days for the IC cooling model to work; otherwise only a fraction 
of shock accelerated electrons will be exposed to IC cooling, and 
the effect described here is weakened. Considering the naturalness 
of this model — it is unavoidable that x-ray photons pass through 
the external shock and cool electrons — and its clear cut, robust, 
implications for external shock lightcurves, one could perhaps turn 
this restriction around and suggest (1) that the angular size of the 
x-ray producing part of the jet increases with time so that it is never 
much smaller than the angular size of the shocked CNM (2) or that 
the external shock does not undergo much sideways spreading for 
at least ~ lO^days even as its Lorentz factor falls with time to order 
unitv^^l. The late time radio lightcurves (fobs ^ 200days) declined 
as ~ i~jj^' (Zauderer et al. 2012; Fig. 3) which is consistent with 
synchrotron radiation from a sub-relativistic blast- wave, and the jet 
spreading is poorly constrained in this case. 

The model predicts a connection between the observed x-ray 



^■^ We don't expect the magnetic field to continue to decline with Lx indef- 
initely. The reason for this is that below a certain threshold x-ray luminos- 
ity magnetic field generation in external forward shock proceeds by Weibel 
mechanism, or some other instability, that is not affected by Lx- 
^^ Kumar & Granot (2003) using 1-D hydro simulations found little evi- 
dence for sideways spreading for a beamed jet interacting with ISM as long 
as it was relativistic. More sophisticated 2-D simulations, e.g. van Eerten 
& MacFadyen (2012), find a similai' result. However, jets appear to un- 
dergo sideways expansion, according to theoretical calculations, when their 
Lorentz factor drops to order unity. 



lightcurve and the temporal and spectral properties of radio and mm 
radiation from Sw J1644-I-57. It requires more extensive numerical 
modeling than presented in this paper to explore this connection in 
detail, which can then be confronted with the extensive data for this 
remarkable event. 

Bamiol Duran & Piran (2013) have suggested two broad class 
of models for Sw JI644-I-57 using equipartition arguments within 
the context of synchrotron emission, but without considering the 
IC cooling of electrons. The first scenario consists of one jet that 
is responsible for both the X-ray and the radio emission. The sec- 
ond scenario invokes two sources: a narrow relativistic jet to pro- 
duce X-rays, and a quasi-spherical, mildly relativistic, source to 
produce the radio emission. Both models require the energy avail- 
able to electrons for synchrotron radiation to increase with time, 
but Bamiol Duran & Piran (2013) argue that the total energy in the 
radio source could be constant. The model suggested in the present 
work roughly corresponds to the first scenario, with a constant total 
energy, where the energy fraction available to electrons for syn- 
chrotron radiation increases with time as the x-ray luminosity de- 
creases and the IC losses weaken. 
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